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A  numerical  technique  has  been  developed  for  capturing  com¬ 
plex,  nonsteady  shock  structures  in  multidimensions.  The 
technique  relies  on  moving  the  computational  mesh  with  the 
shock  wave  so  that  the  features  of  principal  interest  appear 
approximately  stationary.  The  method  has  been  implemented 
using  coordinate-split  Flux-Corrected  Transport  (FCT)  al¬ 
gorithms  which  allow  the  mesh  to  evolve  arbitrarily  with 
respect  to  the  fluid  in  each  coordinate.  The  grid  may  thus 
be  optimized  in  response  to  the  needs  of  a  given  problem. 
Synchronizing  the  grid  and  fluid  motions  permits  significant 
reduction  of  numerical  transients  and  eliminates  numerical 
diffusion.  Shocks  develop  naturally,  with  no  fitting.  The 
method  is  illustrated  by  calculating  complex,  two-dimension¬ 
al  Mach  reflection  phenomena  associated  with  airblasts  and 
shock  diffraction  on  wedges.  The  numerical  results  are  in 
good  agreement  with  available  experimental  data. 

INTRODUCTION 

Numerical  solution  of  transient  multidimensional  gas  dynamics  problems  is 
always  nontrivial.  When,  in  addition,  the  problem  involves  reflecting  super¬ 
sonic  flows,  large  variations  in  length  scales  in  both  space  and  time,  or  phe¬ 
nomena  for  which  neither  analytic  solutions  nor  detailed  experimental  observa¬ 
tions  are  at  hand,  the  state  of  the  computational  art  is  challenged.  Such  a 
problem  arises  in  calculating  the  oblique  reflection  of  shocks  from  solid  sur¬ 
faces  in  planar  geometries  (e.g.  shock  tube  experiments)  or  axisymmetric  geo¬ 
metries  (e.g.  airblasts).  The  complications  arise  mainly  from  the  presence 
of  Mach  reflections  which  occur  when  a  shock  front  impinges  on  a  reflecting 
surface  at  angles  of  incidence  sufficiently  far  from  normal.  The  formation  of 
a  Mach  stem  and,  consequently,  of  a  slip  surface  intersecting  the  triple  point 
(the  confluence  of  the  incident,  Mach,  and  reflected  waves)  results  from  the 
requirement  that  the  flow  behind  the  reflected  shock  be  parallel  to  the  re¬ 
flecting  surface,  which  cannot  be  achieved  through  regular  reflection. 

Attempts  to  calculate  the  properties  of  the  flow  in  Mach  reflections  date 
back  at  least  to  von  Neumann^-  and  the  research  which  grew  out  of  the  wartime 
explosive  studies^-^.  For  the  simplest  problem,  that  of  a  planar  shock 
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reflecting  from  a  plane  surface,  Jones,  Martin,  and  Thornhill  noted  that  it 
is  possible  to  reduce  the  number  of  independent  variables  to  two  by  transform¬ 
ing  to  the  similarity  variables  x/t,  y/t,  a  device  that  was  also  used  by 
Kutler,  et  al®.  Ben-Dor?  developed  a  theory  which  used  shock,  polars  to  explain 
some  of  the  features  of  this  problem,  and  solved  the  system  of  algebraic  equa¬ 
tions  obtained  by  combining  the  jump  conditions  across  the  various  disconti¬ 
nuities  (Courant  and  Friedrichs)®  to  describe  the  flow  in  the  neighborhood  of 
the  triple  point.  To  date,  no  satisfactory  treatment  of  the  complete  flow 
field  has  been  published,  although  some  features  (like  the  shape  of  various 
waveforms)  are  quite  easy  to  model. 

In  connection  with  studies  of  both  chemical  and  nuclear  explosions  there 
have  been  many  attempts  to  model  a  spherical  blast  wave  reflecting  from  the 
ground,  the  so-called  height-of-burst  (HOB)  problem.  The  hydrodynamic  pheno¬ 
mena  in  the  two  cases  are  identical,  although  nonideal  effects  (primarily  ex¬ 
plosive  afterburn  in  the  first  instance  and  radiation  preheating  in  the  second) 
are  different.  Previous  attempts  to  model  two-dimensional  complex  shock  re¬ 
flection  have  suffered  from  restriction  to  describing  part  of  the  system,  the 
use  of  a  special  assumption  like  that  of  self-similarity,  or  less  than  satis¬ 
factory  agreement  with  experimental  data. 9 

The  calculations  discussed  here  represent  a  step  forward  in  overcoming 
these  difficulties.  They  differ  from  previous  numerical  work  in  incorporating 
two  important  computational  developments:  Flux -Corrected  Transport  (FCT)*®  and 
an  adaptive  regridding  procedure,  called  "sliding  rezone", H  which  optimizes 
the  mesh  point  distribution  and  hence  the  resolution  of  surfaces  of  disconti¬ 
nuity. 

FCT  is  a  finite-difference  technique  for  solving  the  fluid  equations  in 
problems  where  sharp  discontinuities  arise  (e.g.  shocks,  slip  surfaces  and 
contact  surfaces).  It  modifies  the  linear  properties  of  a  second-  (or  higher) 
order  algorithm  by  adding  a  diffusion  term  during  convective  transport,  and 
then  subtracting  it  out  "almost  everywhere"  in  the  antidiffusion  phase  of  each 
time  step.  The  residual  diffusion  is  just  large  enough  to  prevent  dispersive 
ripples  from  arising  at  the  discontinuity,  thus  ensuring  that  all  conserved 
quantities  remain  positive.  FCT  captures  shocks  accurately  over  a  wide  range 
of  parameters.  No  information  about  the  number  or  nature  of  the  surfaces  of 
discontinuity  need  be  provided  prior  to  initiating  the  calculation. 

The  FCT  routine  used  in  the  present  calculations,  called  JPBFCT  (an  ad¬ 
vanced  version  of  ETBFCT)^2>consists  0f  a  flexible,  general  transport  module 
which  solves  1-D  fluid  equations  in  Cartesian,  cylindrical,  or  spherical  geo¬ 
metry.  It  provides  a  finite  difference  approximation  to  the  conservation  laws 
of  the  general  form: 


0  f  4>dV  =  -f  <t>  (u-u  )  •  dA  +  J  TdA  (1) 

dt  J  6V(t)  J  bA(t)  "8  */OA(t) 

where  4>  represents  the  mass,  momentum,  energy  or  mass  species  in  cell  6V(t), 
ij  and  u  represent  the  fluid  and  grid  velocities,  respectively,  and  T  repre¬ 
sents  tne  pressure/work  terms.  This  formulation  allows  the  grid  to  slide  with 
respect  to  the  fluid  without  introducing  any  additional  numerical  diffusion. 
Thus,  knowing  where  the  features  of  greatest  interest  are  located,  one  can 
concentrate  fine  zones  where  they  will  resolve  these  features  most  effectively 
as  the  system  evolves  (Fig.  1). 

In  the  next  section  we  describe  the  computational  techniques  used  to  solve 
the  wedge  problem  and  present  the  results  of  four  simulations  carried  out  to 
reproduce  experimental  results  of  Ben-Dor  and  Glass. In  Section  III  we  pre¬ 
sent  a  parallel  discussion  for  a  HOB  calculation.  Finally,  in  Section  IV  we 
summarize  our  conclusions. 
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to  the  reflection  point,  with  stretching  being  initiated  after  the  first  Mach 
stem  reached  'V  20  cells  in  length.  The  timestep  was  recalculated  at  every 
cycle  with  a  Courant  number  of  0.5. 

Figure  2  shows  the  pressure  and  density  contours  and  the  velocity  field  for 
cases  a,b,c,d.  The  pertinent  shock  phenomena  can  be  easily  identified:  inci¬ 
dent  shock,  contact  surface,  first  and  second  Mach  stems.  As  shown  in  Fig.  1, 
the  zoning  is  particularly  sparse  except  for  the  region  of  interest.  Adequate 
resolution  of  the  key  surfaces  (contact  and  second  Mach  stem)  is  obtained  with 
5  zones  in  each  direction.  The  accuracy  can  be  evaluated  by  comparing  the  ex¬ 
perimental  density  distributions  along  the  wall  (Fig.  3). 
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Fig,  2  -  Pressure  and  density  contours  aud  f low  velocity  vectors  (in  frame  of 
reflection  point)  for  planar  waves  with  Mach  number  M  reflecting  from  wedges 
with  angle  0  for  (a)  M=2.03,  0=60°;  (b)  M-2.82,  0-20°;  (c)  M=5.29;  0=30°; 

(d)  M=7.03 ;  0=*5O° 
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Fig.  3.  Comparison  of  density  (in  units  of  ambient  density  p  )  for  cases  (a), 
(b),  (c) ,  (d)  of  Fig.  2  vs.  distance  from  corner.  Points  are  measured  values 
reported  in  Ref.  13. 

HEIGHT  OF  BURST  CALCULATIONS 

Next,  we  performed  a  numerical  simulation  of  a  1KT  nuclear  detonation  at 
31.7  m  HOB,  a  case  which  could  be  readily  compared  with  high  explosive  data.  A 
constant  ambient  atmosphere  was  used  with  a  density  of  1.22  x  10“^  g/cm^  and 
pressure  1.01  x  10^  dynes/cm2.  To  relate  the  energy  and  density  to  the  pressure, 
a  real-air  equation  of  state  (EOS)  was  used.  This  table-lookup  EOS  was  derived 
from  theoretical  calculations  by  Gilmore^’ ^  for  equilibrium  properties  of  air 
and  has  been  vectorized  for  the  Advanced  Scientific  Computer^  .  The  internal 
energy  density  used  in  the  call  to  the  EOS  is  found  by  subtracting  kinetic 
from  total  energy;  this  can  be  negative  due  to  truncation  (phase)  errors.  When 
this  occurred,  the  value  of  the  pressure  was  reset  to  zero. 

The  transition  from  regular  reflection  to  douhle  Mach  reflection  occurs 
at  a  ground  range  approximately  equal  to  the  HOB.  The  size  of  the  mesh  should 
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therefore  be  roughly  twice  the  HOB  in  both  directions.  The  upper  boundary  I 

should  be  far  enough  away  from  the  blast  front  to  be  non-interfering.  We  chose  i 

boundaries  of  55  m  for  the  radial  direction  and  103.5  m  for  the  axial  direction.  ! 

The  fine  grid  in  the  radial  direction  contained  140  out  of  200  total  zones,  ! 

each  5  cm  in  length.  The  rightmost  zones  were  80  cm  in  length,  and  a  smoothing  j 

involving  40  zones  was  performed  between  the  regions  to  guarantee  that  the  zone  j 

sizes  varied  slowly.  In  the  axial  direction  the  fine  grid  contained  75  out  of 
150  total  zones,  each  5  cm  in  length.  Beyond  that  region  tl  e  zones  were  geo¬ 
metrically  increased  by  a  factor  of  1.112. 

Placement  of  the  fine  grid  at  the  origin  of  the  mesh  (ground  zero,  the 
point  at  which  reflection  first  occurs)  was  determined  to  be  optimum  for  cap¬ 
turing  peak  pressure  in  the  airblast  wavefront.  Thus,  as  the  expanding  wave 
moves  along  the  ground  surface,  the  fine  grid  is  always  locked  to  it  and  each 
point  along  the  blast  front  encounters  the  same  spatial  gridding  as  it  approaches 
the  ground.  By  treating  each  point  of  the  incident  front  in  the  same  manner, 
we  insure  that  the  calculation  is  internally  consistent  and  that  the  computed 
transition  point  is  accurate  to  within  the  limits  of  the  resolution. 

The  initialization  provides  a  strong  shock  with  approximate  Mach  number 
M=12.  This  speed  and  the  need  for  restart  capability  led  to  the  choice  of  200 
timesteps  as  an  interval  for  the  spatial  display  (snapshots).  The  dump  inter¬ 
val  that  resulted  was  At  n,  0.3  milliseconds  (ms).  These  dumps  were  stored  on 
magnetic  tape  and  post-processed. 

A  fit  to  the  1-D  nuclear  blast  flow  field  (Ref.  17)  was  used  to  initialize 
the  energy  and  mass  density  and  velocity  field  at  3.76  ms.  The  corresponding 
peak  overpressure  was  113  bars.  After  the  1  KT  flow  field  was  laid  down  inside 
a  radius  of  31.6  m,  the  fine-zone  grid  was  activated  to  follow  the  peak  pressure 
as  it  moved  along  the  ground  surface,  modelled  as  a  perfectly  reflecting  bound¬ 
ary.  This  region  comprised  140  zones,  and  a  switch  was  set  to  keep  40  of  these 
zones  ahead  of  the  reflection  point.  Permeable  boundary  conditions  are  used 
on  the  top  and  right  edges  of  the  mesh,  i.e.,  density,  pressure  and  velocity 
are  set  equal  to  ambient  preshock  conditions.  Reflecting  conditions  were 
applied  to  the  left  and  bottom.  The  total  elapsed  physical  time  in  the  2-D 
calculation,  7.6  ms,  required  5600  cycles.  Times  are  referred  to  t=0  at  the 
start  of  the  calculation. 


The  numerical  simulation  begins  just  before  the  shock  first  reflects  from 
the  ground.  Fig.  4a  indicates  the  pressure  and  density  contours  and  velocity 
vectors  at  time  3.18  ms.  In  Fig.  4b  the  reflected  shock  is  shown  moving  upward, 
the  outward  flow  begins  to  stagnate  at  the  ground  (transition).  Fig.  4c,  t=5.99 
ms,  shows  an  enlargement  of  the  shockfront,  and  the  development  of  the  Mach 
stem,  slip  surface  and  second  Mach  stem.  The  angle  of  the  shock  front  with 
respect  to  the  ground  is  increasing  with  time,  so  that  the  effective  wedge 
angle  is  decreasing.  From  Ben-Dor  and  Glass  1°  one  expects  a  transition  to 
double  Mach  stem  to  occur  at  approximately  45°.  The  angle  in  Fig.  4b  is  about 
45°  and  the  shock  front  has  entered  the  transition  phase.  Figure  4d  shows  the 
fully  developed  shock  structure  at  7.79  ms.  Clearly  visible  is  the  second  Mach 
stem  and  a  vortex  region  behind  the  first  Mach  stem.  Toeing  out  of  the  first 
Mach  stem  can  be  also  seen  in  the  contours  of  Fig.  4d  and  occurs  as  the  fluid 
rolls  forward  where  the  slip  line  would  otherwise  intersect  the  ground.  The 
velocity  field  in  Fig.  4d  also  shows  this  detail. 


One  should  also  note  the  reflected  shock  properties.  The  reflected  shock 
propagates  rapidly  through  the  high  temperature  fireball,  due  to  the  high  local 
soundspeed.  The  shape  of  this  reflected  wave  is  a  primary  difference  between 
the  HOB  case  and  the  wedge  case^.  The  other  major  difference,  of  course,  is 
the  spherically  expanding  blast  wave  which  decreases  in  strength  approximately 
proportional  to  r“  . 
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Fig.  4.  Pressure,  density,  and  velocity  fields  for  HOB  calculation  (a)  in 
regular  reflection  stage;  (b)  at  transition  to  Mach  reflection;  (c)  shortly 
afterward,  when  second  peak  has  become  larger  than  first;  and  (d)  fully  de¬ 
veloped  (note  toe  at  base  of  first  Mach  stem) . 


Finally  we  consider  the  pressure/distance  relation  for  the  HOB  case.  In 
Fig.  5  we  compare  the  results  of  the  numerical  simulation  with  the  data  of 
Carpenter  and  with  empirical  analysis.  Carpenter's  data  are  based  upon  care¬ 
ful  HOB  experiments  with  8  lb  PBX9404  spheres.  The  empirical  analysis  was 
based  on  a  1  KT  nuclear  free  air  curve  and  HOB  construction  factors.  The  cal¬ 
culated  values  in  the  regular  reflection  regime  are  20%  low  and  may  be  attri¬ 
buted  to  a  combination  of  FCT  clipping,  the  resolution  of  the  grid,  and  in¬ 
accuracies  in  the  initialization  of  the  flow  field.  During  and  after  Mach 
reflection,  the  peaks  remain  low  until  the  Mach  stem  structure  has  grown  large 
enough  to  be  resolved  on  the  mesh.  By  the  time  it  occupies  a  region  of  15  cells 
high  and  35  cells  wide,  the  peak  pressures  are  in  good  agreement  with  the  HE 
data  and  the  empirical  analysis. 
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Fig.  5.  Pressure-range  curves  for  first  and  second  (after  transition  -  denoted 
by  TP  -  to  double  Mach  reflection)  peaks. 

SUMMARY  AND  CONCLUSION 

The  complex  2-D  Mach  reflection  phenomena  associated  with  shock  diffraction 
on  wedges  and  height-of-burst  explosions  have  been  modeled  with  the  FAST2D  com¬ 
puter  code.  Four  wedge  cases — regular,  single,  complex  and  double  Mach  reflec¬ 
tion — have  been  calculated  and  the  results  compared  to  experiments.  A  nuclear 
detonation  (1  KT  at  31.7m  HOB)  was  also  simulated.  The  results  give  insight  in¬ 
to  the  formation  and  subsequent  evolution  of  the  Mach  stem,  the  triple  point 
and  the  contact  discontinuity.  The  transition  from  regular  reflection  to  double 
Mach  reflection  is  predicted.  Excellent  agreement  with  Ben-Dor's  data  is  obtained. 
We  suggest  that  the  first  signal  for  transition  is  the  appearance  of  a  second 
peak  behind  the  shock  front  due  to  stagnation  in  the  flow.  Calculated  first  and 
second  pressure  peaks  versus  distance  in  the  HOB  case  agree  both  with  the  HE 
data  and  analysis  to  within  20%. 

The  use  of  the  adaptive  regridding  procedure,  called  "sliding  rezone", along 
with  the  FCT  algorithm  allows  one  to  accurately  predict  the  nonsteady  shock 
structures  in  two  dimensions  for  diffractions  on  wedges  and  HOB  cases.  Compari¬ 
son  with  data  for  both  wedges  and  HOB  yields  the  best  results  obtained  to  date. 
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OICY  ATTN  IN  c-*  4  AN  3? 


ASSISTANT  CHTEr  rr  STAFF 
STUD  I  CS  S  ANALYSES 
OF^A0  TMFNT  OF  THF  A  l  o  frocc 
WASHINGTON,  DC  ?0??0 

OICY  ATTN  AP/SA^I  (TFCH  !.  T«J ) 


ASSISTANT  SECpr-'ADv  rf  TH r  AF 
PFSFA°CH,  OEVF!  f  LOGISTICS 

'1PDAPTMFMT  OF  THF  A|o  FCSCr 
WASHINGTON,  r>o 

OICY  AtTM  S^AlP/rEP  FPP  STP.AT  S  S»ACF  SY« 


BALLISTIC  MissilF  op  p  ICE  /  ‘IN 
A 1 0  FOPCF  SYSTEMS  rC'V  ANC 


N'0?TCM 

Arg,  cm  V>t,cn 
{  V  I  Ml  lT  r  NAN  ) 

OICY 

ATTN 

N  XH 

G  K ALANS* Y 

OICY 

ATTN 

mvmvm 

M  CFivecc.Hin 

OICY 

ATTN 

V  V  V  U 

CFABTPEE 

OICY 

OICY 

ATTN 

ATTN 

vrp'yp 

M' 'V  Y 

0  GAGF 

nF^UTY  rHICF  ctaff 


°FSFA»r.Ht  rFVFl  r  PYF  M  ,  r  ACC 
OF r»ABTMCNT  nF  THr  ATP  cCFCr 
WASt'INGTH  ,  OC  ?n^in 

OICY  ATTN  AFmoj  n  ALFVA^OPOW 
0 1 r  Y  ATTN  ArP  pf>M 
OICY  ATTN  AF  n  OG  I 


CFPAOTVENT  Cc  T HE  AIf>  FCPCF 


DEPUTY  CHIEF  PF  STAFF 
LOGISTICS  C  CMGI^FRING 
HF  P  AP TNENT  OF  THE  AID  ITPCC 
WASH!  NGTC*  ,  DC.  20  3?C 
Olf.V  ATTN'  LFFF 


CrvMAN'QEP 

FP^EIGV  TFCHNOL  nCV  CIV  I  SIGH,  AFSC 
WP  I  GHT-°ATTEP  SON  AFP,  Pi-  45', TT 
OlfY  ATTN  Ml  I S  1  IPP APY 


CPMMAHOFR 

P?^F  A I  P  OE VE  LPPPF^T  rcMrr->t  AF^C 
G°  T  EF  I SS  AF3,  Mv  IT/./,  1 

(  CFS  IpFS  NC  C*"»r  T) 

OlfY  ATTN  TSL r 


STpatcgIC  A  IP  C<"i”wAM? 
r,F°APT^FMT  rif  THF  HIn  FP  PC  E 
rcri)TT  AFB,  W  B  ftp  11  ’ 

OICY  ATTN  NPI-STINFP  U««>/PY 

0 1C Y  ATTN  xprs 

OICY  ATTN  I NT  J  *C'/!NNEV 


VELA  SFISMr»tnCTCAl  C^Trc 
?12  '*3NTGrvEPY  stceft 
AL  E YAPD IPA  ,  VA  ??H4 
OICY  ATTN  G  Mi  p  ICU 
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HFPARTvpait  or  E^POY/OOF  CONTRACTORS 


[)roApr»PMT  OF  EVFROY 
AL  OUClIFPmo  nnFPATTCAS  f'ccjcr 
P  ^  HOY  *5400 
ALTUOnnoniic  f  v»  F 7 1  l  P 

O’CY  ATTN  FTTn 


i*)Fd  AP  np  rvror.Y 

’•>\SHTNOTr»  f  nr  ?0E45 

0 !  C  Y  ATTN  C  p  r  £  T 


OEPAPTMcwr  of  CNFPGY 
NEVADA  norp/vTl^'IS  r:FFICc 
P  0  3°X  l Aim 

ias  vcgas,  nv  mi  14 

OICY  ATTN  f< AT  L  £  PECPPOS  7TS  TFCHMICAl.  LIPP.ATY 


I.AWPCNCF  l  Ivr^HPr  AAT  IO  At  l  At) 

o  n  nny  ppp 
LTVCP’-mpE.  CA  ^4cS0 

OIF Y  ATTN  I  "°C  p  DC‘!G 

OICY  ATTN  I  ?rS  J  FE  ABST  (CLASS 

01C.Y  ATTN  L  <90  n  MFPPTS  (CLASS  L-50'.l 

0  1 C v  ATTN  !  -7  J  KAHN 

n  IT  y  4TTm  n  Gt  rA’N 

OICY  ATTN  L  427  P  SCHCC< 

OICY  ATTN  TCCFNICA!  INFO  OF  *T .  LIP<>APY 
IKY  ATTN  l-TCO  T  t?.lTKOVI  CM 


LOS  ALAMOS  NATTOML  SCTFMTIctp  LAP 
RAIL  STATfrN  mon 
7  n  BOX  1462 
LOS  ALAMOS,  MM  P7S45 

(CLASSIFIED  rN(Y  T  ma?L  STATION  50CO) 
OICY  ATTN  !>  KHTr  jK(:p 
OICY  ATTN  C  YELL  F 0 
OICY  ATTN  M.T.  SANFORD 
OICY  ATTN  •>(  (CLASS  rp:>'’QTS 

OICY  ATTN  E.  JONES 


''-r>flRTMF'lT  OF  FNFPGY/COF  CONTRACTORS 


L07EI  ACE  9IOMPOICAL  r 

cN'-< irnNMFNTAi.  ?<< :y  p*c. 

p  n  npy  «;<?<; o 
41  BUCMFDTHF  ,  \1'  P7I15 

OKY  ATTN  o  ,)ONcS  (  JNCl  fMY) 


HA*'  pnfr  NATrr','t  L  £P.C°A  t^d  V 

NiriFAP  oivism-! 

X-in  LAP  v  f  rnpr*  S  01  VlSir*' 

o  p  npx  x 

PA1'  PlPf.P,  TN  ?7fno 

OKV  ATTN  CIVIL  OFF  PFS  BBOJ 
OKY  Attm  rcflTo.\l  o<",r.H  I.  fP^AFV 


S.\Mrjl\  L A'5n'?ATr'PTCc 
l  T'/cpvpRE  t  ARnt''Ti',,Y 
n  p  f  °  X  Of0 
jri/cp  tppp,  C4 


0 1C  Y  4TTfl  1  IPPA»Y  f  SCC«.*BIT''  CLASS!  rIC  AT  If>  niv. 


SA'niA  NATIONAL  I  AB 
0  0  POX  SPOO 
Al  DUOUEPOStr  »  N*‘  P713S 

(ALL  C’ASS  ATTN  pec  CONTPOt  rFC  FOP  I 
01 CY  ATTN  A  CHA’1A»: 

0 1C  Y  ATTN  |  HU  I 

OKY  ATTN  C°C  I?*"  W  BBPV.F 

01CY  ATTN  A  OAT  M 

OKY  ATT N  M  Prii^PTY 

0 1C Y  ATTN  0141 

01CY  Attn  L  VFPTNAN 

OlCY  ATTN  J.  BANISTER 


f THP°  GOVERNMENT 


CCK>Tt>AL  INTELLir-rvr.r  AGENCY 
WASHINGTON,  nc  70  EG  5 
OICY  ATTN)  r*SWP/*'Fr 


DEPARTMENT  OF  TH r-  lMrP!F7 
RU*FAM  CF  MINES 

01.  nr.  20,  DENVF0  Ff  1EpAl  CENTF* 
DENVER*  CO  00  ??? 

(  (  UNC.L  rvt  V  >  > 

OICY  ATT'.’  TECF  '  IP  {tjNCI  CNIY  J 


DT  p  pC  TOP 

rp-yeg/VL  FMFRC.FNfv  “■  AKAGEPENT  AGENCY 

national  sfc  p*c  -it  if. \?  ion  r.  »sch 

1775  !  ST^EFT, 

WASHINGTON*  DC  2''47'5 

UAU  C  L 4  S c  ATTN  MC*  DOC  CCNTPCL  >=0)1 
OICY  ATTN  •4TT*g*T!CM  K  RGC't  DIV 
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nr  PAP.T^FHT  CF  OfFENSr  CCNT&ACTP^S 


A r  1 1 15 r y  c^p. 

4*5  CLVDr  A VFM‘ !r 
MC: *NT A  r V  v  I  FW ,  CA  0404  > 
Oir Y  ATT\  C  Mfl r 


AfTCS ° ACE  C0°9. 

P  0  ROY  9,t-,57 
L0C.  ANGFl  F<"  t  CA  9000Q 
Olf.Y  ATJN  m  fMPCl? 

01rY  ATTN  TECHNICAL  T  VFr"  VA  T I  ON  S  F  *>  V I C  F  S 


AGP  API  AN  ASSnrTATFS 
N  NASH  ST"  FFT 
FI  SFG'lvon,  C  A  or>->4 5 
O  1 C  Y  ATTN  AGPAPIAN 


analytic  srpvic.cs,  r\c. 

400  APMY-NA Vv  r5P'F 
Amwr.TOf-,  V.A  ??70? 

0  I C  Y  ATTN  0  KCSSFteACHCO 


Aor>L  I  FO  oPSEAOfH  ASSOCIATES,  JVC 
•>401  WYOMING  PLVn  NF  CfTF  j.  _  J 
alrijouf^omf,  m’«  r 7 n r 

01CY  ATTN  J  rPATjcr 
0 1  CY  ATTN  N  HIGGINS 


AO7LTF0  THFn?vt  jnc. 

1010  VFSTWnno  »LVT 
LOS  ANGFI.ES,  CA 

(P  CYS  IF  t  Nf  L  A S S  op  j  cv  IF  Cl  ASS) 
OICY  A TTN  J  TOiiiyr 


■fri  i~y  r  ^ 


PfOAOTMPMT  rc  OEFFNSO  CONTRACTORS 


A0  T  EC  ASSOCIATrc, 

’6046  OOFN  LATINO  or/*o 
MAVWAOO,  CA  <*4*5^5 
''ICY  ATTM  P  MIL 


4STR0N  PESARCH  L  FNGir.FEFIVG 
lOO )  01  n  MIDOLFrinc  K  AY  «1S 
“OJNTAIN  VIEW,  C  A  <5  404  3 
Oir.Y  ATTN  J  Ht'NTTNC-TCN 


r'CV  pE$r ARCH  f.  SYSTEMS  CP.fMp 
’01  LOWFlt  STREET 
WILMINGTON,  MA  01B17 

OICY  ATT*'  1. 1  Hp Ap  Y  AF10 


B0“  COpP. 

7015  JONES  BRANCH  nPI'<E 


“CLEAN, 

VA 

?  ’102 

OICY 

ATTN 

A  IAVACMNC 

OICY 

ATTN 

T  NEJOHPrPF 

OICY 

ATTN 

COR  PC0 ATF  1  f  Cp Apv 

B0M  rnpo, 

R  0  BOX  9.’74 
AL  of O' JEROME  »  NM  F7  11P 
OICY  ATTN  9  HFNSLFY 


rof  t*'G  Co. 

P  0  FOX  1707 
SEATTLE,  WA  P»I?4 
OICY  ATTN  S  STPACK 
OICY  ATTN  APRCCOACF  LIP’APV 
OICY  ATTN  M/S  4? C*1  P  CA* |  SON 
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r»fP4PTvFVT  OF  OEFFNSF  CCNTP&C1PRS 


CM  lrnPMA  prSc.'.r,fH  f  v pr mvp[ toy,  INC. 
626  9  ''i°T  CL  iVf'M  f 
ynooi  AMn  HILL  9,  CA  91? 6 7 
OlfY  4TTN  LIB«'A,’v 
0  If  Y  AITM  K  YRFYELHAGFN 
OlfY  ATTN  m  PC’SFNPIATT 


CALIFORNIA  OFf.F/f'fH  f.  TFf.H*jri  ^9Y,  INC. 
4040  F  I  E  S  T  ST'?  Fc7 

L  JiirpMDPC,  0  4  f'A*’’  Of 

OlfY  ATTN  1  POomal 


CAL  c.  °  AN  f  fo 
n  Q  -IPX  4°0 
B  ’Fr A L C,  NV  l*7?*? 

0  ir y  ATTN  LIf’c,A"V 


CE-NVF&t  UN!Vr°.SITY  rj: 

Colorado  f f m i s p v 

OFNVFP  °FSFAPriJ  INFTITIJT 
»  P  PHY  101  ?7 
OF'NVFP,  CP  30.710 

(ONLY  1  CrPv  CLASS  PpTS) 

01  f.v  AtTN  SFC  nFclCeP  FPp  J  W I  SOTS* I 


coff,  '/.’ASH.  4V  41  v  T  I  f  At  SVC.S  CTP,  INC 

o  f|  pny  10?  lp 
ALPMO’lF7nl»Ef  N’v  9  7  1  I A 


rF**AOTMFNT  OF  OFFFNSC  CCNTPACTCPS 


F»IC  M.  >!ANO 

CIVIL  FV^IMpFPfNjf;  fAC 

UNIVERSITY  OF  mtv  >'Fyjrr 
UNIVERSITY  ST  AT  I rw 
n  q  n  ny 

AL  RtJO'tEPOU!! .  N '*  R  7  1  "* ?. 

0  1C  Y  ATTN  J  l  AN0 
OICV  ATTN  P  LfO-JF 
OICY  ATTN  V  RA"M 
0  ICY  ATJN  J  KQVAPN/ 


GA”0,  IA'C. 

7449  N  NATCHEZ  AVCM  F 
NIL c$  «  IL  60*40 

OICV  AT TM  G  N F l r H A f‘ T  C'<4’C1  ONLY) 


OFNFOAL  EIFCTOTC  CO. 

S".ACF  DIVISION 
VALLEY  FOPGE  Sf*ACr  CFMCP 
n  n  0OX  0556 
OM I L ADEL pH I  A ,  n A  19101 
0  1C. Y  ATTN  M  3CTTNFP 


OFNFRAL  f,CSEA=,('w  C.HPP. 
SANTA  P  A9Q ARA  OTVJ<*,  IfN 
P  0  Rpy  f»7') 

SANTA  PAP HAP A,  CA  93111 
OICY  ATTN  Tfr 


•1-TFCH  LARS,  INC. 

p  n  rrx  1696 

SANTA  MCNTCA,  CA  90406 
OICY  ATTN  3  HAPTEAPAI  v 
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1  i**  ^ «<*■  t-.iv  i  41 


OfTOAQTwpfiT  CF  OFPrnsr  r.fMTRACTOPS 


H'TUZONS  TE r.HMf'i|  r/~,v  t  INC. 
7P30  CLaiPP'^r'iT  Vfsa  PI. VC 

sw  nipr-n,  ca  "mi 

^ir.Y  ATTN  H  KFI'r.c? 


I  I  T  RESEARCH  INSTTTlTr 
1*)  m  ^5TH  ST’>ECT 
CH T  CAG°«  TL  AOA!f 
CITY  ATTN  P  ’  FLO 
01CY  ATTN  '*  JCMNSO 
0  If  Y  AT T N  ry*f  |  im r K‘  T <■,  i_lp&APv 


INF  rt>  MAT  IHN  SCIENCE,  INC. 
123  W  TA3RF  ST*EFT 
SANTA  RAP.  BAP  A,  CA  ‘=3505 
MCY  ATTN  W  n(P7  jjk 


INSTITUTE  rnr  OO^'SF  ANALYSES 
'.00  AP^Y-VAVY  Dp  I V 
A°  L I NOTOV  ,  VA  ?22CH 

O 1 C Y  ATTN  CLASS Tf  IF"  LT3PAPY 


J  0  HALT  TWA  WOO  rCN5t.LT  *NG  SVCS 
p*1  lOAA  CIVIL  CNCTNPFPUn  *1  CO 
N  ppm  IMF  sTr  rcT 
"OA^'A,  ii  ft  ip  01 
OICY  ATTN  s  *1411 


.1.  H.  HI  GO  I  NS 
IASO  S  PACIFIC 
PEOCNOP!  3EACH, 
CICY  ATTN  J 


rr„  INC. 

FCAST  FICFNAY 
CA  FC?77 
CCU  T*S 


4 
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PPPAOTVFMT  Op  OFFFMSf  CCNTRACTf'DS 


K AM AM  4VI0YWF 
9?  SECPvO  AVF'JMP 
mostmwfst  inoustpim  park 

BIPLINOTO*  M  A  PI  -'Pi 


OICY 

ATTN 

P 

P»  c  Tf  \  I  K 

OICY 

ATTN 

l  Tr.PAr  Y 

OICV 

ATTN 

N 

M  r»  r«  *"»  c 

OICY 

ATTN 

r 

C  p  7  s  c  I  r  N 

K  AM  AN  SCIENCES  CC?*1  . 

D  0  nnx  7463 

rrLn!7  4nr'  sppinpc.  cc  anqoo 

OICV  ATTN  0  SACHS 
OtCY  ATTN  F  SHELTCN 
0  7 C  Y  ATTN  L T  Bp  AR  Y 


KAN  AN  TFMPCS 

916  STATE  ST°  EEt  (p  H  C°AWFP  00) 
SANTA  BAPPAP.A,  CA  <7?1P? 

OICY  ATTN  HAS  I  AT 


LDCKHFEP  ’MSS  ILFS  S  SpACF  CC.,  IK. 
o  n  npx  5°4 
SUNNYVAI.Ft  CA  OAF* 6 
OICY  ATTN  .1  WFfSNFP 
0  If  Y  ATTN  TIC  -l  I  GPARY 


MARTIN  M  4°  l  FT  TA  TIP 15  . 
P  0  BOX  5037 
opi.A«nn»  Fi.  37"rs 
OICY  ATTN  0  TCTIFC 


Y.A^TIN  WAP  I  FT  T  4  CC°F. 

o  0  snx  i?<9 

OFNVFP,  CO  on?*' 7 

OICY  4TTN  0  rnpvEP 
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PH0 A°  TMFMT  OF  DFFFNSF  CCNTPACT°RS 


'•nn-Nri.i  pchglas  ccc». 

S  VI  l  B^LSA  AVEM'  e 
H  IVT I NGTr*l  BEACH,  r* 

0  1  C  Y  ATTN  N  HFPrvyv 
0  1C  Y  A T T N  f>  HALJP  IN 
OICY  ATTN  0  PFA*' 


MCOOVMCLL  nn"0lAS  Crco. 
1B5G  I  AKEWCrn  nn;ji_  cy  »c  p 
I  0\G  Br ACH ,  C A  o-3Mf 
0  1  C Y  ATTN  v  pfT'ro 


M  P  E  n  I  T  T  CASES,  !  *'C  . 

?  n  <;ox  i?o.f 
r  e  «>1  AMOS,  CA  V>?7‘» 
OICY  ATTM  J  »nn  jtt 
0  If  V  ATTN  LTPPVJV 


VpTFrP(M,CGY  PCCEA-CH,  T \r. 
46<*  W  wrCnBU°Y 
A1.TA0FMA,  CA  <»  1 0  0  ! 

OICV  ATTN  W  G,,r  r  N 


MISSION  PESFAPfM  CrF’’. 

P  0  OPAWF1'  71 P 

SANTA  BARBARA,  C  A  "  ?  in  ? 

(ALL  CLAPS!  ATTM;  c;  rQ  riFf  FOP  J 
OICY  ATTN  r.  LC?.VMrr 
0 1C v  ATM  0  M rr  *  e  Trc 


°AC IF  If  - S  T rPP  A  °  FSE  AnCH  rrpn. 

ias fj  ci  rvcPFin  t  °i.vc 

SANTA  MCMIC A,  CA  90ACA 
OICY  ATTN  h  ?:PP")E 
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4  v-  feu  4n«w^'.  «.*'»  -x^i-^^w->.i.e  ^na.«iii  >  >1  fcdtMtJENdMlfr  JKtJKb- 


r6«»A0Twt:NT  OF  OFFENSE  CCNTR  ACTORS 


PACIFIC- SIERPA  DESEARCH  roon. 
WAS*<I  NGTO\<  nnpp.STir-NS 
14"!  1  -ML SON  PL''D 
S'JT  Tr  non 

API  T  NGTCN,  V;\  ?  ”* ? 0 *) 

0  1C  v  \TTN  "1  r.rP‘J«LFV 


pact  Fir  a  technology 
P  n  POX  I4R 
DEL  f'AP,  CA 

0  If  Y  ATTN  R  «  JORK 
OICY  ATTN  G  KENT 
Oiry  ATTN  TfPf-  {  fPRAPV 


P  A  T  f-L  cMjrppe  !SrS,  I  AT  . 
p  o  ncv  K31 
HUYTSVILIF,  AL  35E1C 
OKY  ATTN  «  ° A T F L 


PHYSICS  INTFPNAT  IPNAl  CC. 
?' On  ’-*ePCFn  STnrn 


N  LEANDRO, 

CA 

04  f.  7  7 

OK  Y 

ATTN 

1. 

prhrmann 

OICY 

ATTN 

T  r' 

rt-''ICAL  imf>.APY 

OICY 

ATTN 

c 

*Cf’9  r 

OICY 

ATTN 

J 

Tpf",'  sc* 

OICY 

ATTN 

F 

**  AtlFT 

PSD  ASSP'C  I ATCS 
o  n  pox  p6rs 

N  AR  IN  A  CFL  Prv,  r.  A  crpoj 
0 1C Y  ATJN  p  ppRt 
0  1C Y  ATTN  A  KHH! 

OICY  ATTN  J  LFV.^S 
Air.Y  ATTN  V,  ■•’0|-j-T 
OICY  ATTN  J  C  A°  n  F VT  r p 

D1CY  ATT  n  rrrp  f(([  TNF"r>  Y  AT  I  CN  CFNTER 
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nroA*>TMr\T  or  oepcnse  CCNToacto-js 


PANO  COPP. 

170°  vAT'l  STPr"T 
SAk:  TA  mpvjica,  r  A  C04C( 
0  1C Y  ATTN  C  "r'.J 


c.CIEMCr  A^PLIC'C  Tf»**S,  jnc 
3  AO  !ATIGN  I  NS  *J  ye  M  A  t  I  f  f  OTV 
4615  HAWKINS,  Nc 
ALR'JOUFPOUE,  M^irc, 

OICY  ATTN  J  0»S'-T\ 


SCIENCE  APPLICATIONS,  INC. 
p  P  BOY  2351 
LA  JOLLA,  CA  5203° 

OICY  ATTN  w  -IILSCN 

OICY  ATT>'  TFCHICAl  L  !3°  APY 

01CV  ATTN  0  SC»’LAt:C 


sr.irricr  A00l  Ir  AT  TON S ,  INC . 
101  CONTINENTAL  PI  VO 
FI  SEGMNOO,  C  A  PQ74  5 
o ICY  ATTN  0  WOVE 


SCIENCE  A^°LIC.AT  HNS,  INC. 
2450  WASHINGTON  AVEMF 
SAN  lE.V-'OPf',  CA  04577 
OICY  ATTN  3  RFPMSTF  IN 
3 1C  Y  A  T  T  v  3  M  A  y E  l  L 


SCIENCE 
P  0  pox 
MC  L  F  A  N , 

OICY 

OICY 

OICY 

OICY 

OICY 


ADOL I C  AT IC N S ,  INC. 

l?OT 

7  A  ?■*  I  02 

ATTN  J  C OCKAYNP 

AT T*l  B  CHA33FPS  II! 

A TT N  M  KNASFl 

ATTN  W  LAYSCN 
ATTN  P  SIFVEPS 


ppoAPTMCVf  OF  0  c  F  F  M  S  F  CrUTFAFTroS 


Sr!tfwwrST  PrEPAPCI  IN'STiT'JTP 

o  o  D  C  A  V  E  P  '>?ct0 

SAM  ANTCNIH,  TV  7 Pv** 

OICY  ATTN  A  lifNZF). 

O  i  (■  y  4  T  T  N  W  3  A  K  r  S 


SP  I  I  Ni-rrpuftTTC'i'l 
333  F  AVcNSfc'C,r  -A  AVF'MJF 
VEMLC  PAPK«  C A  0  40  ?  5 

OIFY  ATTM  0  A  PP  A  H  A'*  SC  N 


OICY  ATTM  II?PA"V 
OICY  ATTN  J  CCL’CN 


SYSTFYS,  SCIENCE  r  c^rTSAPF  INC 
i>  0  OPT  8243 
AI  BlJ0'IE°0"F  mo  °?1  op 
OICY  ATTN  C  NEE  ^ 


SYSTEMS,  SCIENCE  e  s r F  T w A p F f 
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